245 


SH 2.1-12 


CALCULATIONS OF HEAVY ION CHARGE STATE DISTRIBUTIONS FOR 
NONEQUILIBRIUM CONDITIONS 

A. Luhn, D. Hovestadt 

Max- Planck- 1 nstitut fur extraterrestrische Physik 
8046 Garching, F R G 

ABSTRACT 

Numerical calculations of the charge state distributions of test ions in a hot 
plasma under nonequilibrium conditions are presented In particular, we 
derive the mean ionic charges of heavy Ions for finite residence times in an 
Instantaneously heated plasma and for a non- Maxwell Ian electron 
distribution function. Comparing the results with measurements of the 
charge states of solar energetic particles, we find that neither of the two 
simple cases considered can explain the observations. 

1 Introduction Recent measurements of the charge states of heavy ions in solar energetic 
particle events' indicate that nonequilibrium conditions at the source of the solar energetic 
particles (presumably an active region in the lower corona) are likeley to play a key role in 
establishing the charge state distribution observed at earth. In particular, the temperature as 
derived from the equilibrium tables of charge state distributions 2>3>4 differed significantly 
for different ion species. Calculations of equilibration times for impurity ions m a fusion 
plasma 5 indicate that, ss a function of plasma temperature, the time for initially neutral atoms 
to reach their equilibrium charge state distribution can vary by approximately two orders of 
magnitude We have therefore calculated the mean ionic charges for elements of relevance to 
solar flare studies which are expected a) for a finite residence time of the ions after an 
instantaneous heating of the plasma, and b) for a non-Maxwellian distribution function for the 
electrons m the plasma. 

2 Methods The rate of change of the abundances n, of the ionization states l of an ion with 
nuclear charge Z is given by the following system of equations 

fij = n e (-niSi + n2«2), 

rij = n e ( n,_j S,_j + nj + j «j+j - n,(S, + a,) ), ( 1=2 .Z), (1) 

n z+ i = n e (n z 5 z - n z+l a z+J ). 

For vanishing time dependence this becomes 

nj S, = n l+ | oc , + 1 , ( l = 1 Z) (2) 

Sj denotes the total ionization rate coefficient from ionization state i to l+l, <xj the total 
recombination rate coefficient from i to i- 1 , and n e is the electron density The rate coefficients 

are obtained by multiplication of the relevant cross sections with the electron velocity, averaged 
over the electron distribution function For the time dependent calculations with 8 Maxwellian 
electron distribution we us8 the fit parameters for the rate coefficients given in 2 8nd 3 These 
are calculated using the corona aproximation, i e. neglecting any density dependent terms in the 
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cross sections and photoionization by the radiation field, which Is valid for the parameter range 
considered here The processes included are collisions! ionization and radiative and dielectronic 
recombination together with corrections taking into account autoionization after collisional 
excitation and dielectronic recombination. For the calculation of the charge states produced by a 
non- Maxwellian distribution, however, we need the cross sections a as a function of electron 
impact velocity To obtain these, we use the analytical forms of the rate coefficients for a 
Maxwellian distribution given by 6 and solve numerically for the parameters of the cross 

sections by comparing them with the rate coefficients given by 2 * 3 We then can find the rate 
coefficients for an arbitrary electron distribution f e by calculating Jv8(v)f e (v)dv , were v 

denotes the electron velocity 

The equations ( 1 ) are solved by numerical integration with an algorithm suited for the 
widely differing time scales of the abundance changes for the different ionization states The 
initial charge distributions and the equilibrium distributions for a non-Maxwellisn electron 
distribution function are calculated using the equilibrium equations ( 2). In view of the accuracy 
of the parameters given by 2 for the interpolation formulae for the rate coefficients, conclusions 
based on the values of specific ionization ratios have to be made with caution However, the mean 
values of the charge distributions, which are used for this work, are much less affected 


3 Results and Discussion We first discuss the time development of the mean ionic charge of 
different ions after an instantaneous heating of the plasma. This is motivated by the observation 
that in the initial phase of a solar flare the temperature rises by a factor of approx. 1 0 over the 
ambient coronal temperature of about 10 6 K within a time of the order of 1 second (see for 
instance 7 and 8 ) The subsequent acceleration of protons, alpha particles, and heavier ions may 

remove the ions In a 
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Fig 1 T ime development of mean charges after an instan- 
taneous heating of the plasma from to 6 K to !0 7 K 
at time t=0 


timescale comparable to 
the equilibration time for 
the ionic charge In Fig. 1 

we plot the mean charge 
8s a function of n 0 t for 

the case of a temperature 
rise from 10 6 Kto I0 7 K 
Note that the time for 
Neon and Magnesium to 
reach their new equili- 
brium charge state dis- 
tribution is a factor of 
ten larger than that for 
eg Oxygen if one as- 
sumes a typical electron 
density of 10 9 cm' 3 , it 
is of the order of 1 00 s 
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This is a consequence of the high ionization potential for the K shell electrons of these elements, 
which in this temperature range cannot be removed for ions with a higher atomic number Z The 
differences in the equilibration times contrast with the observation of rather high charge states 
for Ne and Mg 1 , which on the basts of the notion sketched above would be expected for elements 
with short equilibration times. For different initial and final temperatures the equilibration 
times change However, their ratios for the different Ions remain approximately the same We can 
therefore conclude, that the simple model of the freezing in of nonequilibrium charge states due to 
the escape of the particles from the heated coronal material cannot explain the Ionic charges of 
solar energetic particles. 

The second case considered is the calculation of equilibrium charge distributions for a non- 
Maxwellian electron distribution function of the coronal plasma The existence of a non-thermal 
tail of electrons in the transition zone and the lower corona is supported both by theoretical 
calculations for the quiet sun (e.g 9 - 10 ) and observations in flaring regions (eg 1 , -' 2 ) Owocki 
and Scudder ( 6 ) have considered this case and derived the ionization ratios 0 +6 /0 +7 and 

Fe + 1 1 /Fe + 1 2 They find that an estimation of the temperature based on these ratios overestimates 
or underestimates the true temperature of the plasma depending on the ratio kg^/l , were kg is 
the Boltzmann constant, T m the temperature inferred with a Maxwellian distribution, and I the 

ionization potential of the state considered We follow 6 in employing the "kappa" (generalized 
Lorentzian) distribution as a convenient one parameter representation of a Maxwellian with a 
power law tail Our result is that the mean charge changes in general by considerably less than 
the ionization state ratios, especially if states are compared which constitute the flanks of the 
charge distribution 



TelK] 

F ig 2 Mean equilibrium charges of 0 and Mg for a k electron distribution 
function with k= 5 and 10 and a Maxwellian (k=<*>) 

In Fig 2 we compare the mean equilibrium charges of 0 and Mg as a function of temperature 
(defined for the k distribution as corresponding to the mean electron energy) for values of the 
parameter k of 5, 1 0, and (corresponding to a strong and moderately pronounced electron tail 
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and a purely Maxwellian distribution). Fig 3 shows the mean charges for various heavy ions for a 
k distribution with k=5. Although the effect of the energetic tail of the electron distribution is to 
raise the mean charge above that of a Maxwellian with the same temperature, the variation of this 
effect for the different ions is not pronounced enough to explain the variation of the temperature 
inferred for the different ions from a Maxwellian 



Fig 3 Mean eguihbnum cnarges as a function of electron tempereture for 
a k distribution with k=5 

In closing we conclude that neither of the above two simple models is able to explain the 
observed discrepancy in the temperatures derived from the mean ionic charge states for the heavy 
constituents of flare generated solar energetic particles. 
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